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ABSTRACT 

The clustering properties of moderate luminosity {Lx = 10^^ — 10^^ erg s^^) X-ray 
selected AGN at 2 « 0.1 are explored. X-ray sources in the redshift interval 0.03 < z < 
0.2 are selected from a serendipitous XMM survey of the SDSS footprint (XMM/SDSS) 
and arc cross-correlated with the SDSS Main galaxy sample. The inferred X-ray AGN 
auto-correlation function is described by a power law with amplitude vq « 5h^^Mpc 
and slope 7 « 2.0. The corresponding mass of the dark matter haloes that host X- 
ray AGN at 2: « 0.1 is « 10^^ Mq. Comparison with studies at higher redshift 
shows that this mass scale is characteristic of moderate luminosity X-ray AGN out 
to z w 1. Splitting the AGN sample by rest-frame colour shows that X-ray sources in 
red hosts are more clustered than those associated with blue galaxies, in agreement 
with results at 2 « 1. We also find that the host galaxies of X-ray AGN have lower 
stellar masses compared to the typical central galaxy of a w 10^^ Mq dark matter 
halo. AGN hosts either have experienced less stellar mass growth compared to the 
average central galaxy of a w 10^^ Mq halo or a fraction of them are associated 
with satellite galaxies. 

Key words: galaxies: active, galaxies: haloes, galaxies: Seyfert, X-rays: diffuse back- 
ground 



1 INTRODUCTION 

Observations (e.g. Ferrarese & Merritt 2000, Gebhardt et 
al. 2000) and theoretical arguments (e.g. Silk 1998, King 
2005) suggest that the growth of supermassive black holes 
(SMBHs) at the centres of galaxies are intimately related 
to the formation of their hosts. However, the details of the 
interplay between these two components are still not well 
understood. 

The environment of AGN, i.e. the mass of the typical 
dark matter halo (DMH) in which they live, is a powerful 

diagnostic of the physics that drive the formation of SMBHs 
and their hosts. Different models for the co-evolution of 
AGN and galaxies make specific testable predictions on how 
the environment of active SMBHs depends on redshift and 
accretion luminosity (e.g. Marulli et al. 2006, Hopkins et al. 
2007, Fanidakis 2010, Degraf et al. 2010, Bonoh et al. 2009). 
Moreover, there has been progress recently in phcnomeno- 
logical methods to associate baryons to DMHs at different 
redshifts, to infer in a least model dependent way how galax- 
ies evolve at different environments (e.g. Zlieng et al. 2007, 
Conroy & Wechsler 2009, Zehavi et al. 2011, Avila-Rccsc 
& Firmani 2011). By determining the typical DMH of AGN 
one cam place them in the context of those studies to explore 
the conditions under which SMBH grow at different epochs. 



An advantage of studying the environment of AGN is that 
it provides one of the few diagnostics of the AGN/galaxy in- 
terplay that is immune to contamination of the stellar light 
of galaxies by emission from the central engine. AGN can 
easily outshine their hosts, thereby rendering studies of the 
stellar mass, star-formation history and morphology of their 
hosts challenging and prone to systematics. This poses a se- 
rious limitation in the comparison between observations and 
models for the growth of SMBHs and highlights the impor- 
tance of clustering studies. 

Powerful UV bright QSOs are one of the few AGN 
classes for which tight constraints on their large scale distri- 
bution are available. This is because of the apparent bright- 
ness of these sources in the optical which has allowed large 
spectroscopic follow-up programs, such as the 2QZ (2dF 
QSO Redshift Survey, Groom et al. 2004), SDSS (Sloan Dig- 
ital Sky Survey, Richards et al. 2002, Schneider et al. 2005) 
and 2SLAQ (2dF-SDSS LRG and QSO Survey Cannon et 
al. 2006, Groom et al. 2008). The picture emerging from 
those studies is that powerful UV bright QSOs live in dark 
matter haloes of few times 10^^ Mq almost independent 
of redshift and accretion luminosity (e.g. Groom et al. 2005, 
Myers et al. 2007, da ingela et al. 2008, Ross et al. 2009, 
Ivashchenko et al. 2010, but see Mountichas et al. 2009). 
These properties are broadly consistent with the predictions 
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of the merger model for the growth of SMBH (e.g. Hopkins et 
al. 2007, Bonoh et al. 2009), suggesting that powerful QSOs 
are the product of interactions between gas rich galaxies. 

The QSOs selected in the above surveys however, repre- 
sent luminous and rare sources above the knee of the AGN 
luminosity function. Therefore, they are not representative 
of the overall AGN population and their contribution to the 
accretion power of the Universe is small. Additionally, their 
selection at UV/optical wavelengths raises concerns on pos- 
sible biases against obscured sources. Observations at X-ray 
wavelengths provide an efficient way of selecting AGN over a 
wide luminosity baseline nearly independent of obscuration. 
Study of the clustering of X-ray AGN can therefore con- 
strain the fueling mechanism of the sources that dominate 
the accretion history of the Universe. Despite considerable 
progress in the last few years however, the large scale dis- 
tribution of X-ray AGN is still not well constrained. This is 
primarily because of the apparent faintness of these sources 
at optical wavelengths and their relatively high surface den- 
sity on the sky, which make large spectroscopic follow-up 
programs, similar to those carried out for UV bright QSOs, 
extremely expensive in telescope time. 

In the absence of any follow-up observations, the most 
widely used approach for studying the clustering of X-ray 
AGN is the angular auto-correlation function (e.g. Basilakos 
et al. 2004, 2005, Plionis et al. 2008, Ebrero et al. 2009). The 
main limitation of this approach is that assumptions have 
to be made on the redshift distribution of the AGN popu- 
lation to infer the mass of their dark matter haloes. This 
introduces systematics and model dependent biases which 
can be largo. This problem has been mitigated by extensive 
follow-up spectroscopic programs in selected few X-ray sur- 
vey fields (e.g. Barger et al. 2003, Brusa et al. 2009), which 
allowed estimation of the real-space auto-correlation func- 
tion of X-ray AGN samples (e.g. MuUis et al. 2004, Gilli 
et al. 2005, 2009, Yang et al. 2006, AUevato et al. 2011). 
One of the results from these studies is the importance of 
sample variance. Spikes in the redshift distribution of X- 
ray sources can significantly affect clustering studies in cur- 
rent X-ray surveys, which typically have small angular sizes 
(e.g. Gilli et al. 2005, 2009). An alternative approach for 
studying the large scale distribution of AGN is the real- 
space cross-correlation function with galaxies over the same 
cosmological volume (e.g. Coil et al. 2009, Krumpe et al. 
2010). This approach has certain merits compared to auto- 
correlation function methods, as long as the clustering of 
the galaxy population is known to a good level of accu- 
racy and spectroscopy is available for the AGN. Firstly, the 
space density of galaxies is typically much larger than that 
of AGN, thereby suppressing random errors when counting 
AGN/galax;y pairs. Secondly, by cross-correlating AGN with 
galaxies in the same field the impact of sample variance is 
minimised. The studies above suggest that X-ray AGN live 
in dark matter haloes with masses 5x 10^^ — 5x 10^^ h~^ M©, 
which are, on average, more massive than those of UV bright 
QSOs. The wide range in the estimated DMH masses could 
be the result of random errors and systematics that affect in- 
dividual measurements or because the AGN clustering may 
depend on the accretion luminosity, the redshift and/or the 
level of obscuration of the central engine (e.g. Plionis 2008, 
Hickox et al. 2011, AUevato et al. 2011). 

In this paper we explore changes in the large scale dis- 



tribution of moderate luminosity {Lx < 10*^ ergs ^) X-ray 
selected AGN from z ~ 1 to z « 0.1. We use a serendip- 
itous XMM survey of the SDSS area (XMM/SDDS, Geor- 
gakakis & Nandra 2011) to compile a sample of low red- 
shift AGN, z ~ 0.1. These sources are then cross-correlated 
with the SDSS Main Galaxy sample (Strauss et al. 2002) 
to infer their clustering properties. The advantage of the 
serendipitous XMM/SDSS survey is that the AGN selection 
function at low redshift is almost identical to that of X-ray 
AGN at z ~ 1 detected in deep Chandra and XMM surveys. 
Therefore differential selection effects between low and high 
redshift samples are minimal. This allows direct compari- 
son of the environment of moderate luminosity X-ray AGN 
across redshift to investigate possible evolutionary trends. 
Throughout this paper we adopt Ho — 100 kms~^ Mpc~^, 
= 0.3 and = 0.7. Rest frame quantities (e.g. lu- 
minosities, dark matter halo masses) are parametrised by 
h = Ho/wo. 



2 THE X-RAY AGN SAMPLE 

The clustering properties of X-ray AGN at z ~ 0.1 
are investigated by selecting low redshift sources detected 
in a serendipitous XMM survey of the SDSS footprint 
(XMM/SDSS). The construction of the XMM/SDSS source 
catalogue, including X-ray source detection, flux estimation 
and optical identification, is described in Georgakakis & 
Nandra (2011). 

The XMM/SDSS survey includes pointings targeting 
clusters of galaxies. The ovcrdensity of sources in those 
fields may bias large scale structure studies. Therefore 
XMM/SDSS survey fields that have clusters as their prime 
targets are excluded from the analysis. These observations 
are identified from the target name keyword of the event 
files. This reduces the total XMM/SDSS survey area to 
102 deg^ 

The low redshift X-ray subsample of the XMM/SDSS 
survey consists of 175 serendipitous hard-band (2-8 keV) and 
297 full-band (0.5-8 keV) detections with 0.03 < z < 0.2, 

X-ray luminosity Lx(2 — lOkeV) > 10'*^ erg s^^ (see bo- 
low) and r < 17.77 mag after correcting for Galactic extinc- 
tion (Schlegel et al. 1998). The magnitude cut corresponds 
to the limit of the SDSS Main Galaxy Sample (Strauss et 
al. 2002), which provides the majority of redshifts in the 
SDSS. The photometry is from the New York University 
Value-Added Galaxy Catalog (NYU-VAGC, Blanton et al. 
2005) which corresponds to the SDSS DR7 (Abazajianet al. 

2009) . The X-ray luminosity cut is to limit the sample to 
AGN which contribute substantially to the X-ray luminos- 
ity density of the Universe at low redshift (e.g. Aird et al. 

2010) . As discussed by Georgakakis et al. (2011) contami- 
nation by normal galaxies is not a concern for luminosities 
Lx(2-10keV) > lO-^^crgs-i. 

Clustering results arc presented for both the full and 
the hard-band selected AGN samples. The former has the 
advantage of larger size, thereby improving the statistical 
reliability of the results. The latter is selected at rest-frame 
energies of about 2-9 keV at 2 = 0.1, which are similar to 
those of Chandra and XMM X-ray AGN samples at « « 1 
(typically l-14keV). This facilitates the comparison across 
redshift by minimising differential selection effects. 
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X-ray luminosities are estimated in the 2-10 keV band 
after correcting the observed flux in the 2-8 keV band (hard- 
band selected sample) or the 0.5-8 keV band (full-band se- 
lected sample) for intrinsic absorption parametrised by the 
hydrogen column density, Nh- F'or individual X-ray AGN 
this quantity is determined from the hardness ratios between 
the soft (0.5-2 keV) and the hard (2-8 keV) X-ray bands as- 
suming an intrinsic power-law X-ray spectrum with index 
r = 1.9 (Nandra et al. 1994). 



3 METHODOLOGY 

3.1 AGN clustering estimation 

Because the low redshift subset (0.03 < z < 0.2 and 
r < 17.77 mag) of the XMM/SDSS survey source catalogue 
is small, we choose to quantify their clustering properties by 
estimating their cross-correlation function with the much 
larger sample of SDSS spectroscopically identified galax- 
ies. We use the NYU-VAGC to select a total of 592,017 
sources in the SDSS Main Galaxy spectroscopic sample (bit- 
mask parameter vag_SELECT equals 7) with redshifts in 
the interval 0.03 — 0.2 and extinction corrected magnitudes 
r < 17.77 mag. 

Incompleteness in the galaxy redshift catalogue because 
of fiber collisions is a source of bias in clustering studies. 
The SDSS fibers have finite size and cannot be placed closer 
than 55arcsec, which corresponds to about 70/i~'^kpc at 
z = 0.1. This separation is too small to affect our results 
and conclusions as the X-ray/galaxy cross-correlation signal 
is dominated by pairs on much larger scales (see Results 
section). We nevertheless correct for this effect by assigning 
a source that was not observed because of fiber collisions the 
redshift of the galaxy with which it collided, as proposed by 
Blanton et al. (2005). 

The determination of the cross-correlation function 
with a large sample of galaxies is superior to the estima- 
tion of the auto-correlation function in the case of small 
samples, like X-ray AGN, because random errors are sig- 
nificantly suppressed. An additional advantage is that the 
cross-correlation requires knowledge of the selection func- 
tion of galaxies only, which is less complex than that of X- 
ray AGN. Also, sample variance is affecting in the same way 
X-ray AGN and galaxies. The impact of this bias is therefore 
minimised in the cross— correlation function calculation. The 
cross-correlation approach however, requires an accurate es- 
timate of the auto-correlation function of galaxies to infer 
the clustering properties of AGN. 

Next we present the equations used to determine the 
clustering of AGN. Since both the auto-correlation and the 
cross-correlation functions are special cases of the 2-point 
statistics of the AGN and galaxy populations, they are both 
defined by the same basic equations. In this section the term 
correlation function refers to either the auto-correlation or 
the cross-correlation functions. When necessary we will dif- 
ferentiate between the two quantities. 

The real space correlation function, ^, is calculated as 

Nrd DD 



Ngal DR 



1, 



(1) 



r, i.e. AGN-galaxy pairs in the cross-correlation function 
or galaxy-galaxy pairs in the case of the auto-correlation 
function, DR are the AGN-random pairs (cross-correlation) 
or galaxy-random pairs (galaxy auto-correlation function) 
at separation r. Random points within the surveyed area 
are produced by randominsing the positions of SDSS galax- 
ies taking into account the SDSS window function and the 
spectroscopic completeness of the galaxy sample. The ran- 
dom catalogues provided as part of the NYU-VAGC data 
release version 7.2 are used. They contain ~ 3x more ran- 
doms than galaxies. The random points are distributed with 
constant surface density in the SDSS footprint as defined by 
the Large Scale Structure mask of NYU-VAGC (Blanton et 
al. 2005). The redshifts assigned to the randoms points fol- 
low the redshift distribution of the galaxy sample. 

When the correlation function is measured in redshift- 
space, the clustering is affected at small scales by the rms 
velocity dispersion of AGN along the line of sight and by 
dynamical infall of matter into higher density regions. If si 
and S2 are the distances of two objects 1, 2, measured in 
redshift-space, and 6 the angular separation between them, 
then a and n are defined as 



TV = (s2 — Si), along the line-of-sight, 
= ^ ^ ^ — —S, across the line-of-sight. 



(2) 
(3) 



These are small angle approximations. Equation [T] then be- 
comes: 

Nrd DD{g,-k) 

? tt) = — — — — - 1, (4) 

Ngal DR{a,Tx) 

To the first order, the non-linear redshift-space distortions, 
i.e. the small-scale peculiar velocities, appear only in the 
radial component. These effects are therefore minimised by 
integrating ^ along the vr direction. The resulting two-point 
statistic is the projected correlation function 



^(cr, n)dTV. 



(5) 



where Nrd is the number of random points, Ngai is the num- 
ber of galaxies, DD are the data-data pairs at separation 



In practice the maximum scale of the integration is set to 
TVmax = 70h~^Mpc (da yingela et al. 2008). If larger scales 
are included, the signal will be dominated by noise. If the 
integration is limited on small scales the amplitude will be 
underestimated. 

In the estimation of the projected AGN/galaxy cross- 
correlation function each AGN is weighted by the factor 
1/Vmax, the inverse of the maximum volume within which 
a source can be detected. This accounts for the complex 
selection function of the XMM/SDSS AGN. Because of vi- 
gnetting and the spatially varying width of the XMM's Point 
Spread Function (PSF) the solid angle within which the sur- 
veyed area is sensitive to X-ray faint sources is smaller than 
for brighter ones. As a result lower luminosity AGN are un- 
derepresented in the low redshift subset of the XMM/SDSS 
survey, but in a way that can be accurately quantified and 
corrected for in the Vmax calculation. It is emphasised that 
this is a second order correction and does not have a strong 
impact on the results and conclusions. For simplicity this 
weight is not included in the equations presented in this 
section. 

Under the assumption that the real-space correlation 
function follows a power-law of the form ^(r) = 
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the real-space correlation length, ro, and slope, 7, can be 
estimated directly from the projected correlation function 



Wp{a) 



r(|)r( 



7 — 1 > 



r(i) 



(6) 



where r(a;) is the Gamma function. This is because equation 
[5] can be rewritten as 



Wp((T) 



rC(r) 



idr, 



(7) 



and then solved analytically by substituting the power-law 
form for ^(r). 

The projected correlation function can also be used to 
estimate the real-space correlation function, ^(r), even if a 
power law form is not adopted. By inverting Wp{a) (Saun- 
ders et al. 1992) 



1 



duj{a)/da 



da 



and assuming a step function for uip(a) = Wi it is found 



1 \ ^ LU 

n ^-^ a 



-ln{- 



-) 



(8) 



(9) 



for r = Oi. That measurement is generally more noisy than 
Wpio). The projected correlation function and equation [6] 
are therefore used to infer ro and 7, under the assumption of 
power-law for 5. However, we check that the ro and 7 values 
obtained by equation |6] provide a good approximation to the 
real-space correlation function infered from equation [9] 

The above procedure allows us to measure the AGN- 
galaxy projected cross-correlation function, Wp[AG), and 
the galaxy projected auto-correlation function, WpiGG). Us- 
ing these measurements and assuming a linear bias, the AGN 
projected auto-correlation function, Wp{AA), is 



Wp^AA) 



Wp{AGf 



(10) 



Wp{GG) 

For the estimation of errors of the correlation function mea- 
surements the survey area is split into six subregions, each 
of which includes nearly equal number of X-ray AGN. The 
cross-correlation function is then estimated for each of the 
six subregions. The variance at a given scale is 

JV 

(11) 



iV _ 1 ^ DR 

L = l 

where N is the number of fields, i.e. N=6, DRi^ is the AGN- 
random pairs in the field, DR is the overall number of AGN- 
random pairs, is the correlation function measured in a 
subregion and ^ is the overall correlation function. 

3.2 Estimation of the bias parameter 

The bias parameter of AGN relative to the underlying dark 
matter halo distribution is estimated using two different ap- 
proaches that are often adopted in the literature. This is 
to facilitate direct comparison of our results with previous 
studies. 

The first method uses the integrated projected 
AGN/galaxy cross-correlation function and normalizes the 
result to the volume contained in a sphere with radius of 



20/1"^ Mpc (e.g. Ross et al. 
da Angela et al. 2008) 

3 
20^ 



^(r)r dr. 



2009, Mountrichas et al. 2009, 



(12) 



where ^20 is the integrated cross-correlation function. The 
lower and upper limits are set to 1 and 2Qh^^ Mpc, respec- 
tively. For the estimation of ^20, our uip(a) measurements 
are used. The galaxy bias, bo, is estimated as 



1,2 



£,GG 



bo 



fzGG 
?20 



en 
S2 



(13) 



where ^20™ is the integrated correlation function of dark 
matter. It is estimated from the normalized ACDM power 
spectrum model of Smith et al. (2003) for cosmological pa- 
rameters n„^{z = 0) = 0.3, q.k(z = 0) = 0.7, r = 0.17 

and (Ts = 0.8, in accordance with the recent WMAP results 
(Spergel et al. 2007). The value of erg has been revised by the 
latest analysis of the WMAP data from 0.84 to 0.8. To facil- 
itate the comparison with previous studies we also present 
results for as = 0.84. 

Having measured the galaxy bias we can then estimate 
the AGN bias, bA, from the relation 



bAba 



^AG 



bA 



CAG 
S20 



bG C2" 



(14) 



where ^20 is the integrated AGN/galaxy cross-correlation 
function. The equation above also assumes scale indepen- 
dent bias. 

The second approach for estimating the bias uses the 
best-fit ro and 7 values of the inferred AGN auto-correlation 
function, obtained from equation 10 (e.g. Krumpe et al. 
2010). The clustering strength is expressed in terms of the 
rms fluctuation of the density distribution over a sphere with 
a comoving radius of 8h~^Mpc 



'^S.AGiV = "^2(7) 



where 



^0 



8/i-iA/pc 



72 



"^'^^^ (3-7)(4-7)(6-7)2^' 
the AGN bias can be calculated via 



bAGI^ 



as, AGN 



(15) 



(16) 



(17) 



The bias of AGN or galaxies is related to the mass of the 
dark matter halos they live in (e.g. Mo & White 1996). 
Therefore their clustering properties can be used to infer 
their dark matter halo masses. In this calculation the el- 
lipsoidal collapse model of Sheth, Mo & Tormen (2001) is 
adotped. The methodology described by da Angela et al. 
(2008) and van den Bosch (2002) is followed to convert the 
bias measurements to dark matter halo masses. We use the 
bias values estimated from the as,AGN rms fluctuation of 
the density distribution. 



4 RESULTS 

Figure [l] (left panel) compares the full-band AGN-galaxy 
cross-correlation function, Wp{AG), with the galaxy auto- 
correlation function, Wp{GG). Galaxies are less clustered 
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Figure 1. Left panel: The full-band AGN-galaxy cross-correlation is shown by filled circles. Open circles show the galaxy auto-correlation 
results. Solid triangles connected with the solid line show the inferred full-band AGN auto-correlation function. Right panel: Same as in 
the left panel but using the hard-band AGN sample. For clarity, open circles and solid triangles are offset in the horizontal direction by 
(5 logo- = —0.02 and -1-0.02, respectively. 



than AGN at all scales. The infered AGN projected auto- 



correlation function, Wp{AA), (i.e. equation 10 1 is also plot- 
ted in the figure. The errors in Wp{AA) are estimated 
by adding in quadrature the uncertainties of Wp{AG) and 
■Wp{GG). The right panel of Figure [l] shows the clustering 
results for the hard-band AGN sample. Within the errors the 
clustering of those sources is similar to the full-band sam- 
ple. This is also shown in Tables [l] and [2] which present the 
best-fit estimates for the ro and 7 for the cross-correlation 
and the auto-correlation functions respectively, assuming a 
power-law form for ^. Scales in the interval 0.25 — 17h~'^Mpc 
are used to fit the data. The errors correspond to the 68th 
percentile around the minimum . 

The assumption that the real-space correlation function 
is well approximated by a power law is justified by Figure 
[3] It compares the inferred ^(r) for the AGN/galaxy cross- 
correlation and the galaxy auto-correlation function without 
making any assumptions on its functional form (i.e. equa- 
tions Is] ^ with power-law fits using the ro and 7 values of 
Tables[iland[2l 

We can therefore use the best-fit parameters for ro and 7 
to estimate the bias of AGN and galaxies (equations 12 - 17 1. 
The results are presented in Table|3j where h{^2o) and b(cr8) 
are respectively the biases estimated from the first and the 
second method described in section 3.2. The errors are deter- 



mined from the variance across the six subregions of the sur- 
veyed area. Within the uncertainties the bias values, calcu- 
lated from the two different methods are consistent. We cau- 
tion that the h{(Ts) and b(^2a) are determined from slightly 
different scales. The former is based on the ro and 7 pa- 
rameters estimated by fitting the projected cross-correlation 
function on scales of 0.25 — 17h~^Mpc. The latter is deter- 
mined by integrating the ^(r) from 1.0 to 20h~^Mpc. We 
have confirmed however, that this difference does not change 
the results and conclusions. 

Next, we explore variations in the clustering of X-ray 
AGN with rest-frame optical colour. This is motivated by 
previous studies suggesting that AGN in red cloud hosts are 
more clustered than those in blue cloud galaxies (e.g. Coil 
et al. 2009). For this exercise we use the full-band AGN 
sample because of its larger size. Red/blue galaxies/AGN 
are separated by their "'^(u — g) rest-frame colour, i.e. the 
difference between the absolute magnitudes of the source in 



the" 



0.1 



g bands, which are the SDSS u, g filters shifted to 



z = 0.1. The calculation of ^'^{u — g) is carried out using the 
KCORRECT version 4.2 routines (Blanton & Roweis 2007). 
The distribution of X-ray AGN and galaxies in "'^(it — g) 
colour is plotted in Figure [4] The division between red and 
blue galaxies and AGN is set to "^'^(u — g) ~ 1.5 (Blanton 
2006). There are 325,510 galaxies and 154 X-ray AGN with 
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Table 1. rg and 7 values for f(r) (scales of 0.25 — 17/i~^Mpc) for the cross-correlation measurements. 



AGN (full ) 


Red AGN 


Blue AGN 


AGN (Hard) 


Red galaxies 


Blue galaxies 


4.4+0- 
7 2.00iro^ 


4.9 ±0.4 
1.98 ±0.10 


4.1 ± 0.2 
1.95 ±0.05 


4.3 ±0.2 
1.90 ±0.11 


4.80 ±0.06 
1.82 ± 0.01 


3.60 ±0.05 
1.71 ± 0.01 



Table 2. ro and 7 values for ^(r) (scales of 0.25 — 17h ^Mpc) for the galaxy auto-correlation and for the inferred auto-correlation 
function using equation |10| 





AGN (full ) 


Red AGN 


Blue AGN 


AGN (Hard) 


galaxies 


Red galaxies 


Blue galaxies 




5.0 ±0.5 


6.6 ±0.6 


4.7 ±0.4 


4.8 ±0.6 


4.00 ± 0.05 


5.00 ±0.10 


2.90 ± 0.05 


7 


2.00 ±0.11 


1.82 ± 0.19 


2.12 ± 0.16 


2.02 ±0.18 


1.76 ±0.01 


1.90 ±0.02 


1.58 ±0.01 



• Full-band AGN autocor 




1 10 
CT(h"' Mpc) 

Figure 2. Filled circles show the inferred full-band AGN auto- 
correlation and open circles the inferred hard-band AGN auto- 
correlation. For clarity, open and filled circles are offset in the 
horizontal direction by Slogcr = —0.01 and +0.01, respectively. 



colours bluer than that cut (blue subsamples) and 266,508 
galaxies and 143 X-ray AGN with °-'^{u - g) > 1.5 (red 
subsamples). 

Figure [5] plots the projected cross-correlation of red and 
blue galaxies with the full galaxy sample. The relevant ro 
and 7 best-fit values, bias parameters and dark matter halo 
masses are shown in Tables[2]and[3] As expected red galaxies 
are more clustered than blue ones at all scales, in agreement 
with previous studies (e.g. Madgwick et al. 2003, Zehavi et 
al. 2005, Coil et al. 2008, Hickox et al. 2009). The projected 
cross-correlation function of red/blue AGN with the overall 
galaxy population is shown in Figure |6] The clustering re- 
sults are presented in Tables [2] and [S] X-ray AGN follow the 



same pattern with galaxies, i.e. active SMBH in red hosts 
are more clustered than those in blue galaxies. 

This is further demonstrated in Table |4] which presents 
the relative bias of different subsamples, defined as the 
ratio between the Wp{a) measurements of the subsample, 
fere! = \/ 'Wp{a)i / Wp{a)2 (Coil et al. 2007), integrated over 
two scales, 0.25 - 8h"^Mpc and 1 - 8h"^Mpc. The errors 
on these measurements are estimated from the variance of 
the relative bias across the six subregions of the survey area 
(Section 3.1). This table confirms that red AGN and galax- 
ies cluster more than blue AGN and galaxies, respectively. 
Also, the overall X-ray AGN population has similar clus- 
tering properties as red galaxies but is more clustered than 
blue ones. 

We also investigate the dependence of the clustering 
on luminosity and X-ray obscuration by splitting the sam- 
ple into two nearly equal size groups at logLx = 41.8 
(erg/s) and logA'^^f = 22 (cm~^) respectively. We do not 
find any statistically significant trends of the AGN clus- 
tering with luminosity or X-ray hardness. This null result 
maybe be because of the relatively narrow luminosity and 
obscuration baselines of our low redshift X-ray AGN sam- 
ple. Cappelluti et al. (2010) for example, find a higher clus- 
tering for type I AGN compared to type lis. Their sam- 
ple is drawn from the SWIFT-BAT AGN catalogue, which 
is sensitive to heavily obscured systems. Also, Krumpe et 
al. (2010) find evidence for higher clustering only for very 
bright AGN, Lx > 10** ergs "\ The low redshift subset of 
the XMM/SDSS survey does not have the volume to detect 
such powerful sources. 

Finally, we calculate the masses of the dark matter halos 
that host AGN and galaxies, as described in Section 3.2. The 
results appear in Table [3] We do not estimate a lower limit 
for the mass of the halos that host blue galaxies. This is 
because the Icr lower limit of the bias of those sources is 
fe = 0.69 and the minimum bias value derived from Sheth et 
al. (2001) is fe = 0.72. The calculations show that moderate 
luminosity AGN at 2 ~ 0.1 reside in dark matter halos with 
typical mass Mdmh « lO" h"^ Mq. 

The AGN bias and DMH mass estimates presented in 
Table 3 are robust to the details of the adopted methodol- 
ogy. Cosmic variance is a concern as in the cross-correlation 
function estimation we use all SDSS galaxies, not just those 
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Figure 3. Left panel: The full-band AGN-galaxy cross-correlation function in real-space is shown by filled circles. The solid line shows 
the fit on scales 0.25 < r < f7h^^Mpc. Open circles show the galaxy auto-correlation results and the fit is shown by the dashed line. 
Right panel: Same as in the left panel but using the hard-band AGN sample. For clarity, open circles are offset in the horizontal direction 
by 5 log cr = -1-0.03. 



overlapping with the XMM pointings. Therefore the cosmo- 
logical volumes of the AGN and galaxy samples are different. 
We repeat the analysis for the full-band selected AGN sam- 
ple using only those SDSS galaxies that lie within the XMM 
fields of our serendipitous survey. This excerise yields the 
same AGN bias as in Table 3. Additionally, the exclusion 
of XMM pointings targeting clusters (see section 2) might 
appear arbitrary and subjective. We include AGN in those 
fields and cross-correlate them with the SDSS galaxies in 
the same cosmological volume, i.e. those overlaping with the 
XMM pointings of the XMM/SDSS survey. The infered full- 
band AGN bias is 10 per cent lower compared to that listed 
in Table 3, i.e. within the estimated errors. 



5 COMPARISON WITH PREVIOUS STUDIES 

In this section our bias estimates for X-ray AGN at z « 0.1 
are compared with previous studies that select active SMBH 
either at X-rays or UV/optical. To facilitate the direct com- 
parison and present the results from different papers in a 
uniform manner we re-estimate the bias from each study us- 
ing the quoted ro and 7 values based on the ag methodology 
of section 3.1. Studies that have measured only the redshift- 
space correlation function and do not provide real-space ro 
and 7 values are excluded from the analysis. Moreover, an- 



gular auto-correlation function studies are not included in 
our compilation as the deprojection of the angular signal to 
3-dimensions using the Limber's formula (Limber 1953) may 
introduce systematic uncertainties. 

Table [S] presents the studies used for these measure- 
ments and the most important properties of the samples 
they used. The errors on b(ag) are calculated based on the 
errors on ro, 7. The bias parameter is compared with our 
measurement in Figure [7] and the expected evolution of the 
bias for different dark matter halo masses. Below we com- 
ment on selected datapoints appearing in Table |5] 

From the high redshift samples plotted in Figure [7| the 
Coil et al. (2009) datapoint is the most appropriate to com- 
pare with our measurement of the bias aX z ~ 0.1. They 
select X-ray AGN at 2; ~ 1 at rest-frame energies of about 
1-14 keV similar to those used here. The luminosity range 
of their AGN is similar to our low redshift sample. They 
also use the cross-correlation function with galaxies to im- 
prove on the statistics and minimise sample variance effects 
in their calculations. The mass of the dark matter halos that 
host these high redshift X-ray AGN, as estimated based 
on the quoted ro and 7 values by Coil et al. (2009), is in 
good agreement with our measurements at low redshift. This 
shows that moderate luminosity X-ray AGN live in dark 
matter haloes of similar mass, ~ lO^^ft-^Mg, at all red- 
shifts out z ~ 1. 
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Table 3. Bias values for AGN and galaxy samples using the two approaches discussed in the text. Values marked with * are calculated 
for as = 0.84. The ^20 values are from the cross-correlation measurements (except for galaxies). The M^omh are based on h{as = 0.8). 
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Figure 4. The O i?i— g histograms for the galaxies (solid line) and 
the AGN (dashed line). For clarity, the AGN histogram has been 
scaled, by the ratio Nq/Na, where Nq and Na are the number 
of galaxies and AGN, respectively. 



In agreement with Coil et al. (2009) we also find that 
AGN in red hosts are more clustered than those in blue 
galaxies. However, they find that the difTerence in AGN clus- 
tering due to the host color is more apparent for their 2 ~ 1 
AGN than in our sample at 2 = 0.1. The relative bias be- 
tween red and blue AGN at 2 ~ 1, brei ~ 1-9 (see Table 2 
of Coil et al. 2009) compared to bmi ~ 1-3 at 2 ~ 0.1 (Table 
[4|. Similarly, our relative bias estimation for the red/blue 
galaxies is 6rei ~ 1.35. Coil et al. (2009) do not give a value 
for this measurement, but using their quoted relative bias 
values for red/blue AGN, red AGN/red galaxies and blue 
AGN/blue galaxies we infer a red/blue galaxy relative bias of 
brei = 2.02 ±0.49 at 2 « 1. This apparent discrepancy is be- 
cause of the differential evolution with redshift of the bias of 
dark matter halos of different masses. This is demonstrated 



in Fig.|8]which shows how the relative bias of red/blue AGN 
and galaxies evolves with redshift. 

Gilli et al. (2005) and (2009) measured the projected 
auto-correlation function of X-ray AGN in the Chandra 
Deep Field South (CDFS), Chandra Deep Field North 
(CDFN) and the COSMOS surveys. AGN in the CDFs have 
luminosities comparable to those of Coil et al. (2009), al- 
though the COSMOS sample is more X-ray luminous, be- 
cause of the large angular size (2deg^) and shallower depth 
of the X-ray survey in that field. GiUi et al. (2005, 2009) find 
that sample variance has a strong effect on their results, even 
in the wide-angle COSMOS sample, resulting in large values 
for the clustering length and the corresponding dark matter 
halo masses, logMoMH ~ 13.5 h~^ ^0- Narrow peaks in the 
AGN redshift distribution have a significant contribution to 
the clustering signal and are primarily responsible for the 
high Mdmh estimates. Preferentially removing galaxies in 
those peaks decreases the AGN bias in the CDFS and the 
COSMOS to values comparable to those of Coil et al. (2009) 
and to our estimates. However, MaruUi et al. (2009) argue 
that sample variance cannot account for the large clustering 
length in the CDFS. They use semi-analytic models to follow 
the cosmological evolution of AGN and produce mock AGN 
catalogues with selection functions similar to the CDFN and 
the CDFS. Although the clustering of AGN in their simu- 
lations is in good agreement with the observations in the 
CDFN (Gilli et al. 2005) they cannot account for the large 
clustering signal in the CDFS within 2.0 — 2.5cr. These re- 
sults highlight the importance of developing methods which 
account for sample variance (e.g. cross-correlation function, 
see also AUevato et al. 2011). 

In addition to X-ray selected AGN, Figure [t] (crosses) 
and Table [5] also include measurements for the bias of 
UV/optically selected QSOs from Mountrichas et al. (2009), 
Ross et al. (2009) and Ivashchenko et al. (2010). The latter 
study measured the auto-correlation function of ~ 50, 000 
SDSS quasars at 2 = 1.47. We estimate that their quasar 
bias is & = 2.06 ± 0.10 which is higher than their calcula- 
tion, b = 1.44 ± 0.22. This discrepancy is likely because of 
the different methodology they follow for the bias estima- 
tion. They use the redshift-space and real-space measure- 
ments to calculate the infall parameter (/3; see their eqn 7) 
and based on this calculation they estimate the AGN bias 
(/9 = fi^m, z)/b; see their section 4.1 for more details). Ross 
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Figure 5. Cross-correlation functions of blue galaxies {^'^{u — 
g) < 1.5) with the galaxy sample (filled triangles) and red galax- 
ies with the galaxy sample (open triangles). Red galaxies cluster 
more than blue galaxies. 



et al. (2009) used ~ 30, 000 spectroscopic SDSS quasars, 
over the redshift range 0.3 ^ z ^ 2.2. We estimate a bias for 
their QSOs of & = 2.30i;j;2o (Table [Hjl. This is higher than 
b = 2.06±0.03 that they find using a methodology similar to 



that of Ivashchenko et al. (2010). The discrepancy between 
the two numbers is likely because of the different approaches 
used to infer the QSO bias. 

Figure [7] shows that in a redshift range < 2 < 2 
the AGN bias datapoints cluster around the curve that 
correponds to dark matter halos with mass Mdmh ~ 
lO^^h~^M0. This seems to be the case for both moder- 
ate luminosity X-ray AGN and powerful UV/optical quasars 
as presented above. We caution however, that many optical 
studies (Croom et al. 2005, Myers et al. 2007, da Angela 
et al. 2008) find lower masses for the dark matter halos of 
powerful UV bright QSO, Mdmh ~ 5 x lO^^h'^Mg. The 
results from those studies are not shown in Figure [7] for the 
reasons discussed at the beginning of the section. 



6 COMPARISON WITH MODELS OF AGN 
FUELING AND EVOLUTION 

This paper explores the environment of moderate luminosity 
X-ray AGN at 2 ~ 0.1 by estimating their cross-correlation 
function with the SDSS Main Galaxy spectroscopic sam- 
ple. The adopted methodology has advantages compared to 
auto-correlation function studies, in that random and sys- 
tematic uncertainties can be better controlled. Our cluster- 
ing methodology, AGN selection function and X-ray lumi- 
nosity range are very similar to those of Coil et al. (2009) at 
2 ~ 1 . This allows direct comparison of the clustering prop- 
erties of moderate luminosity AGN across redshift by min- 
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Figure 6. Blue {°-'^{u-g) < 1.5) AGN cross-correlation function 
with galaxies (filled triangles) versus the red-AGN/galaxy cross- 
correlation function (open triangles). Red AGN cluster more than 
blue AGN. For clarity, open and filled triangles are offset in the 
horizontal direction by 5 log a = —0.01 and -1-0.01, respectively. 



imising differential selection effects. Our results combined 
with those of Coil et al. (2009) show that moderate lumi- 
nosity X-ray selected AGN live in dark matter halos with 
masses Mdmh ~ 10^'^ h~^ M© at all redshifts since z « 1. If 
powerful UV/optically selected QSOs reside in lower mass 
halos (few times lO^^h~'^M0), as suggested by some stud- 
ies (e.g. Croom et al. 2005, Myers et al. 2007, da Angela 
et al. 2008, Ross et al. 2009, but see Figure then the 
fueling mode of those sources may be different from that of 
the X-ray selected moderate luminosity AGN studied in this 
paper. 

Semi-analytic models for the growth of SMBH that as- 
sume mergers as the main mechanism for triggering AGN 
activity, predict parent dark matter halos similar to those 
determined for UV/optically selected QSOs and lower than 
those estimated here for moderate luminosity X-ray AGN 
(e.g. Marulh et al. 2008, Bonoh et al. 2009). Thus our re- 
sults argue against major mergers as the main channel for 
fueling the SMBH in moderate luminosity AGN. AUevato et 
al. (2011) have recently expanded this conclusion to higher 
X-ray luminosities, Lx ~ 10** ergs~^. They used the XMM- 
COSMOS field to estimate dark matter halos in excess of 
10^^ h"'^ Mq for powerful X-ray sources. This suggests that 
a substantial fraction of the accretion density in the Uni- 
verse is associated with dark matter haloes with masses 
~ 10^"^ h~^ M0, higher than what merger models predict. 

Stochastic accretion as described in Hopkins & Hern- 
quist (2006) cannot explain the massive dark matter halos 
of moderate luminosity AGN. In this model, disk instabil- 
ities or minor interactions fuel at high accretion rates rel- 
atively small SMBHs in spiral galaxies with abundant cold 
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Figure 7. Comparison of our bias estimation (filled circle) to other studies. Crosses present the estimated biases for optical samples. 
Prom low redshift to high redshift, Mountrichas et al. (2008), Ross et al. (2009) and Ivashchenko et al. (2010). The dotted lines present 
the expected b{z) of dark matter halo masses. Dark matter halos have been estimated as described in the text (in units of h~^Mpc). 



gas supply. This mechanism produces by design moderate 
luminosity AGN, which however, are predicted to lie in low 
density environments, as this is where large gas reservoirs, 
i.e. blue galaxies, are typically found, particularly at low 
redshift. 

A more promising alternative is the model of Ciotti & 
Ostriker (2001), in which stellar winds from evolved stars 
in quiescent early type galaxies provide the fuel to supply 
the SMBH. The accretion energy heats the ambient gas. 



slowing down subsequent infall when the Compton tem- 
perature of the emitted radiation is higher than the mean 
galactic gas temperature. Evolution in this case is character- 
ized by strong oscillations, in which very fast and energetic 
AGN bursts are followed by longer periods during which 
the SMBH is dormant. This fueling mode applies to ellipti- 
cal galaxies, which reside, on average, in dark matter halos 
similar to those of our moderate luminosity AGN sample. 
One approach for getting insights into the growth of 
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Table 4. Results for the relative bias. 



Sample Relative bias Relative bias 
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Figure 8. Relative bias for red/blue AGN (filled circle at z = 0.1) 
and red/blue galaxies (filled triangle). Triangles are offset by 0.01 
for clarity. Dotted lines show how the (relative) bias evolves with 

redshift. At z = 0.9 the br-^i of Coil ct al. (2009) are plotted. The 
higher relative bias values at z = 0.9 are in agreement with those 
at 2 = 0.1, given the redshift evolution of the relative bias. 



SMBH is to investigate how the baryonic matter evolves 
within the typical dark matter haloes that AGN are found. 
Empirical methods have been developed recently to connect 
galaxies to haloes at different epochs to constrain the his- 
tory of galaxy assembly (e.g. Zheng et al. 2007, Conroy & 
Wechsler 2009, Zehavi et al. 2011, Avila-Reese & Firmani 
2011). Although these methods cannot directly probe the 
formation of SMBH at the centres of galcixies, they can pro- 
vide useful information on the conditions under which black 
holes grow. 



A key result from the studies above is that the ratio 
of galaxy to halo mass shows a peak at the dark matter 
mass scale of about 10^^ h"'^ M©, which depends only mildly 
on redshift, and then decreases at higher/lower masses (e.g. 
Leauthaud et al. 2011, Conroy & Wechsler 2009). The halo 
mass where star-formation is most efficient is close to that 
measured for luminous UV/optically selected QSOs, which 
are proposed to be the products of gaseous major mergers 
(e.g. Di Matteo et al. 2005, Hopkins et al. 2008, MaruUi et 
al. 2008, Bonoli ct al. 2009). In contrast, moderate lumi- 
nosity X-ray AGN at both z ^ 0.1 and 2 ~ 1 (e.g. Coil 
et al. 2009) live in haloes that are, on average, offset from 
the sites where star-formation proceeds more efficiently. The 
typical dark matter halo of these systems (lO'^^ h"'^ M©) cor- 
responds to the scale where satellites start to dominate the 
stellar mass budget of the halo (e.g. Leauthaud et al. 2011), 
i.e. small and moderate size groups (e.g. Coil et al. 2004). 
This might indicate that moderate luminosity X-ray AGN 
may be fueled by the cooling of hot gas associated with the 
group environment (e.g. Croton et al. 2006, Bower ct al. 
2006). Semi-analytic models that include prescriptions for 
the growth of SMBHs predict that this accretion mode can 
produce luminosities comparable to those of the AGN sam- 
ple presented here (e.g. Fanidakis ct al. 2011). Alternatively, 
the increasing importance of satellites in the integrated stel- 
lar mass of lO'^^ h~^ M© haloes suggests that they may play 
a role in the activation of the SMBH, e.g. via tidal interac- 
tions or mergers with the central galaxy of the halo. Semi- 
analytic models of galaxy evolution (e.g. Bower et al. 2006, 
De Lucia & Blaizot 2007) predict that the main route of 
stellar mass growth since ^; = 1 of the central galaxies of 
haloes with present-day mass 10^^ h~^ M© are mergers with 
satellites (Zehavi et al. 2011). As long as AGN are associated 
with the central galaxy of a halo (e.g. Starikova et al. 2010), 
such gravitational interactions may also be responsible for 
triggering accretion of material onto the central SMBH. 

We caution however, that studies that connect stellar 
and halo masses in a phenomenological way find contradic- 
tory results on the importance of mergers in haloes with 
present-day mass lO^^h'^M©. Conroy & Weschsler (2009) 
argue that in such haloes star-formation is primarily respon- 
sible for the growth of the stellar mass of the central galaxies 
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Table 5. Clustering measurements for X-ray and UV bright AGN from the literature. Columns are: (1) reference to the AGN sample. 
A star marks UV/optical selected QSOs; (2) methodology used to determine the clustering, i.e. auto-correlation or cross-correlation; (3) 
name of the survey that the AGN sample was selected from; (4) number of sources used: (5) scales used each study to determine the 
clustering amplitude and slope; (6) best-fit clustering amplitude measured in each study; (7) best-fit clustering slope measured in each 
study; (8) mean redshift of the sample; (9) the bias values re-calculated using the crs methodology discussed in the text. 



Study 
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Methodology 
(2) 
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(3) 


no 

of objects 
(4) 
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(h-iMpc) 
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'■{) 
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7 
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z 

(8) 


b(CT8) 
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0.25 


- 17 
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2.00± 0.11 


0.10 
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^■^'-'-0.17 
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cross-cor 


XMM/SDSS 


175 


0.25 


— 17 


4.8 ±0.6 


2.02 ±0.18 


0.10 
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Coil et al. (2009) 


cross-cor 


AEGIS 


113 


0.1 


- 8 


5.95 ±0.90 


1.66 ±0.22 


0.90 


1 q7+0-26 
-0.25 


Krumpe et al. (2010) 


cross-cor 


RASS 


1552 


0.3 


— 15 


4.281^5^ 


^•"'-0.12 


0.27 
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.^•J^J^_0.12 


MuUis et al. (2004) 
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NEP 


219 


5.0 


- 60 




1 85+1-^° 
—0.80 


0.22 




GiUi et al. (2005) 


auto- cor 


CDFS 


97 


0.2 


- 10 


10.3 ± 1.7 


1.33 ±0.14 


0.84 


2.64 ±0.51 


Gilli et al. (2005) 


auto- cor 


CDFN 


160 


0.2 


- 10 


5.5 ±0.6 


1.50 ± 0.12 


0.96 


1.87 ± 0.25 


Gilli et al. (2009) 


auto-cor 


COSMOS 


538 


0.3 


- 40 


8.65i;;:i;, 


i.88_(, (,7 


0.98 




Yang ct al. (2006) 


auto-cor 


CDFN 
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0.2 
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t.3»_o.l4 


0.8 
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Yang et al. (2006) 
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Starikova et al. (2010) 


auto-cor 


Bootes 


1282 


0.5 
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5.4 ±0.5 


1.97 ±0.09 


0.37 


1.52 ±0.16 


Starikova et al. (2010) 


auto-cor 


Bootes 
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0.5 
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7.0 ± 0.8 


1.97 ±0.09 


1.28 


2.98 ± 0.36 


Cappelluti et al. (2010) 


auto-cor 


BAT 


199 


0.2 - 


-200 




1 64+°-°* 


0.045 




Koutoulidis et al. (in prep.) 


auto-cor 


CDFs/ECDFs 
COSMOS/AEGIS 
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0.2 
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6.1 ±0.2 


1.8 (fixed) 


1.24 


2.4 ±0.1 


Mountrichas et al. (2008)* 


cross-cor 


2SLAQ 
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5.0 - 


- 25.0 


6.7 ±0.6 


1.75 ± 0.10 


0.55 


1.67 ±0.29 


Ross et al. (2009)* 


auto-cor 


SDSS 
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1.0 - 


- 25.0 
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1 90+°°* 
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Z.OU_g 20 


Ivashchenko et al. (2010)* 


auto-cor 


SDSS 
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1.0 


- 35 
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1.87 ±0.07 


1.47 


2.60 ±0.17 



since z ~ 1. In contrast, Zheng et al. (2007) find that mergers 
with smaller galaxies, eitiier satellites or central galaxies of 
smaller haloes, drive the stellar mass assembly since z w 1 in 
those haloes, in qualitative agreement with SAMs (Zehavi 
ct al. 2011). The discrepancy is likely related to the way 
galaxies are associated with dark matter haloes in different 
studies and the assumptions adopted to determine the con- 
tributions of star-formation and merging to the stellar mass 
assembly. 

A limitation of schemes that connect galaxies to dark 
matter haloes is that they describe the average stellar mass 
growth as a function of both galaxy and halo mass. There is 
however, evidence that moderate luminosity AGN hosts are 
different from the typical galaxy that lives in haloes sim- 
ilar to those of AGN. Using the results of Behroozi et al. 
(2010, their Table 3) we find that dark matter haloes with 
masses of ~ lO^''h~^M0 at z = 0.1 and z = 1 are pre- 
dicted to have central galaxies with mean stellar masses of 
about 5 X lO^°h"^M0 and 6 x 10^°h"^MQ, respectively. 
The scatter in stellar mass at a fixed halo mass is about 
0.16 dex (Behroozi et al. 2010). Moderate luminosity X-ray 
AGN hosts however, have a broad stellar mass distribution 
(e.g. Shi et al. 2008, Bundy et al. 2008, Georgakakis et al. 
2011). Using the stellar mass function of X-ray AGN from 



Georgakakis et al. (2011) for example, we estimate mean 
stellar masses of about 1.5 X lO^°h"^M0 at a « 0.1 and 
2 X 10^°h"^M© z w 0.8 with a variance of about 0.4 dex 
at both redshifts. This indicates that X-ray AGN, for their 
dark matter haloes, live in smaller than average galaxies in 
terms of stellar mass. X-ray AGN hosts have therefore ex- 
perienced less merging and/or less star-formation (possibly 
due to feedback processes associated with the central en- 
gine) compared to typical galaxies in dark matter haloes of 
similar mass. Alternatively, this might be interpreted as ev- 
idence that a fraction of the X-ray AGN are not associated 
with the central galaxies of their haloes but with satellites, 
which are expected to have a much wider stellar mass func- 
tion (e.g. Leauthaud et al. 2011). Recent studies on the Halo 
Occupation Distribution of AGN suggest that at least some 
of them live in sattelite galaxies (e.g. Padmanabhan et ak. 
2009, Miyaji et al, 2011). Contrary to that, Starikova et al. 
(2010) find that X-ray AGN reside close to the centres of 
their dark matter haloes and tend to avoid satellites. 

Evidence that X-ray AGN are in lower mass galaxies 
for their dark matter haloes already exists in the literature. 
Coil et al. (2009) find that galaxies hosting X-ray AGN are 
more likely to reside in more massive dark matter halos than 
the overall galaxy population of the same color and optical 
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luminosity. Digby North et al. (2011) also found that X- 
ray AGN live on average in higher density environments 
compared to stellar-mass matched galaxy samples. 



7 CONCLUSIONS 

The XMM/SDSS survey is used to constrain the cluster- 
ing of moderate luminosity (Lx(2 — lOkeV) = 10*^ — 
10**ergs~^) X-ray selected AGN at z ~ 0.1. These sources 
are found to reside in haloes with mass of about ~ 
10^^h~^MQ. Comparison with studies at higher redshift, 
suggest that this halo mass scale corresponds to the typical 
environment of X-ray AGN in the luminosity interval above 
at least out to 2 ~ 1. 

Haloes with masses ~ lO^'^h'^M© correspond to the 
group scale, where satellites dominated the stellar mass bud- 
get of the halo. This suggests that either accretion of hot gas 
associated with the group environment or interactions of the 
central galaxy with satellites are responsible for the activa- 
tion of the SMBH in our AGN sample. 

Splitting the AGN sample by colours shows that those 
with red colours are more clustered than those in blue 
colours. This is in agreement with results at higher redshift, 
2 ~ 1. We do not find a dependence of clustering on obscu- 
ration or accretion luminosity. A wider luminosity baseline 
than that spanned by our sample is probably needed to ex- 
plore trends of the AGN clustering with Lx ■ 

There is also evidence that the hosts of moderate lu- 
minosity X-ray AGN have lower stellar masses compared to 
the typical galaxy in haloes of the same size. This may have 
important implications for understanding the conditions un- 
der which supermassive black holes at the centres of galaxies 
grow. 
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